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We use the tight-binding model and the random-phase approximation to investigate the 
intrinsic plasmon in silicene. At finite temperatures, an undamped plasmon is generated 
from the interplay between the intraband and the interband-gap transitions. The extent of 
the plasmon existence range in terms of momentum and temperature, which is dependent 
on the size of single-particle-excitation gap, is further tuned by applying a perpendicular 
electric field. The plasmon becomes damped in the interband-excitation region. A low 
damped zone is created by the field-induced spin split. The field-dependent plasmon spec¬ 
trum shows a strong tunability in plasmon intensity and spectral bandwidth. This could 
make silicene a very suitable candidate for plasmonic applications. 
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I. INTRODUCTION 


A carbon nanosheet, namely graphene, was discovered a few years ago. The unique low- 
lying linear energy dispersions have provided researchers with an abundance of new physics, 
including the Klein paradox for tunneling through a rectangular electrostatic potential 
barrier BEE], the anomalous quantum Hall effect BE], and peculiar optical properties 
!3 Hums Ho]. Following the successes achieved with graphene, researchers have intensified 
their search for other graphene-like two-dimensional materials HHH21H3]. In recent years, 
silicene, a nanosheet consisting of silicon atoms, has been synthesized unmans!. Belonging 
to the same group as graphene in the periodic table, silicene is predicted to exhibit the 
feature-rich electronic properties. Furthermore, it could be better suitable for practical 
applications than graphene for two main reasons. The first is its compatibility with current 
Si-based device technologies. The second is its sizable energy gap, owing to the larger 
magnitude of its atomic spin-orbit interaction (SOI) and its buckled structure. The buckling 
means that the two sublattices are separated by a vertical distance of 21 (« 0.46A), which is 
a result of the larger ionic size of silicon atoms compared to carbon atoms. Other group-IV 
crystals, e.g., Ge and Sn, are also expected to appear in this form in their 2D structures 

mm- 

The out-of-plane buckling suggests that an on-site potential difference between the A 
and B sublattices is tunable under a perpendicular electric held [T9, 20]. This allows for 
the control of the band gap, along with the spin and valley polarized states. Accord¬ 
ingly, many spin and valleytronics, like the quantum spin Hall effect |2T, 221 E3J and the 
optically-generated spin-valley-polarized charge carriers [54] [25], 26j [27], may be investi¬ 
gated. Additionally, when the electric potential exceeds the magnitude of the SOI, 2D 
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silicene is predicted to undergo a topological phase transition from a topological insulator 
to a band insulator, a property that should be important in electronics and optics. 

Collective-Coulomb excitations, dominated by the electron-electron interactions, are 
useful for understanding the behavior of electrons in a material due to the involvement of 
screening effects. This basic property has been widely studied in various low-dimensional 
systems, e.g., carbon nanotubes [2BJI2H], 30], graphene layers [ 311 13211331131115511351 137 ] . and 
other novel materials [38] [39, 30]. In monolayer graphene, plasmons (quantized collective- 
excitation modes) hardly exist at low frequencies due to the lack of free carriers. This may 
be improved by doping or insertion of a gate to increase the free-charge density inra,i.e., 
changing the Fermi level in an extrinsic condition. Alternatively, an intrinsic low-frequency 
plasmon (with zero Fermi energy) may be induced by increasing the thermally excited elec¬ 
trons and holes in the conduction and valence bands, respectively [ 33 ]. This is based on the 
strong dependence of the carrier density on the temperature (n oc T 2 ). The temperature- 
induced plasmon is also predicted to exist in monolayer silicene [31] • The difference is that 
the temperature-induced plasmon in graphene is always located in the interband region and 
suffers the Landau damping [ 33 ]. However, in silicene, the temperature-induced plasmon 
could be in a single-particle-excitation (SPE) gap arising from the spin-orbit interaction 
[43] • The opening of the SPE gap makes the plasmon undamped. In this paper, we investi¬ 
gate the effects of a perpendicular electric field on the T-induced plasmon in silicene. The 
energy-loss function is used to derive the plasmon spectrum. Three plasmon modes with 
distinct degrees of Landau damping are found in different momentum, temperature, and 
field strength ranges. The occurrences of plasmons and their discontinuous dispersions are 
associated with the special structures in the dielectric functions. 
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II. METHODS 


Similar to graphene, silicene consists of a honeycomb lattice of silicon atoms with two 
sublattices made up of A and B sites. The difference is that silicene has a buckled structure, 
with the two sublattice planes separated by a distance of 21 with / = 0.23 A (Fig. 1). In the 
tight-binding approximation, the Hamiltonian for silicene in the presence of SOI is written 


as [?, [35]: 


H 



t , ■^ so 

,ai ° 3V3 


«ij»a/3 
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( 1 ) 


where c\ a (typ) creates (annihilates) an electron with spin polarization a (/ 3 ) at site i ( j ). 
The sum is taken over all pairs of nearest neighbor (i,j) or next-nearest-neighbor (( i,j )) 
lattice sites. The first term in Eq. (1) is the usual nearest-neighbor hopping integral with 
a transfer energy t — 1.6 eV. The second term describes the effective SO coupling with 
Aso=3.9 meV. a = (a x , a y , a z ) is the vector of the Pauli matrices. Vij is determined by the 
orientation of the two nearest bonds connecting the next-nearest neighbors, i.e., = +1 

(—1) if the next-nearest-neighbor hopping is anticlockwise (clockwise) with respect to the 
positive z axis. The third term is the intrinsic Rashba SO coupling with Ar 2 = 0.7 meV, 
where dA = dij/\dij\ with dij connects two sites i and j in the same sublattice, and Uij = ±1 
represents the A(B) site. The fourth term is the staggered sublattice potential energy 
produced by the external electric field, where /i, = ±1 for the A(B) site. The tight-binding 
Hamiltonian can also be used to describe other group-IV elements, such as germanium with 
a different set of parameters (t = 1.3 eV, Aso=43 meV, Ai?2=10.7 meV, and 1=0.33 A in 
the case of Ge USD- 

When the system is perturbed by a time-dependent Coulomb potential, the electron- 
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electron interactions would induce a charge density fluctuation which acts to screen the 
external perturbation. This dielectric screening determines the normal modes of the charge 
density oscillations. Electronic excitations are characterized by the transferred momentum 
q and the excitation frequency oj, which determine the dielectric function. Within the 
random-phase approximation, the dielectric function can be expressed as: 

e(q,u) = e 0 -v q x°(q,uj) ( 2 ) 

, where v q = 2ite 2 /q is the in-plane Fourier transformation of the bare Coulomb potential 
energy, and e 0 =2.4 (taken from graphite [46]) is the background dielectric constant due to 
the deep-energy electronic states. It is noted that changing eo would make a vertical shift of 
the real part of dielectric function (Fig. 3) and thus change the peak intensity and position 
in the energy loss function (Fig. 4). However, the plasmon behavior is qualitatively the 
same (Fig. 5). Such a form of dielectric function is utilized in many researches on two- 
dimensional systems, theoretically (311 [4?, 0SJ and experimentally [4S : , EH, El] • The bare 
response function y° is given by [46, 52]: 


x°(q,u) 
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E 

n,n r 


IstBZ 


dk x dk y 


| (n'\ k + q\e l ^' r \n] k) | 2 


f(E n '(k + q))-f(E n (k)) 
E n ' (k + q) — E n (k ) — {u + zT) 


(3) 


Each Bloch state is labeled by the Bloch wave vector k and the band index n. E n is 
the corresponding eigenvalue. The Fermi-Dirac distribution is f(E n ) = l/{ 1 + exp[(E n — 
/i)//csT]}, where ks is the Boltzmann constant and fi is the chemical potential, p is 
set to zero over the range of the symmetric conduction (c) and valence (u) bands under 
the investigation. Y is the energy broadening which results from various de-excitation 
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mechanisms, e.g., the optical transitions between the valence- and conduction-band states. 
It is usually set small in the low-frequency region. Changing h would not alter the peak 
positions in the energy-loss function (Fig. 4), but the peak intensities and widths will be 
affected. 

The details of the calculation of Coulomb matrix elements are shown below. 

\(n'] k + q\e l ^' r \n] k) = 

Y ' si (k + q)u* nsi {k)\. (4) 

s=a,/3 i=l,2 

T\ and T 2 define the positions of atoms in a unit cell. u ns i(k ) {u n ' si (k+q)) are the coefficients 
for the TB wave functions derived from Eq. (1). {(f> z (r — fi)\e~^'^ r ~ Ti ' > \(f> z (f — f^)) = I(q) = 
[1 + l 2 ^] 2 ]” 3 was calculated by using the hydrogenic wave function [46], where a 0 is the 
Bohr radius and Z is an effective core charge [53]. For small q’s, I(q) is very close to 1. 
The double integral is performed numerically with respect to wave vector and band index. 

III. RESULTS AND DISCUSSION 

The Hamiltonian of Eq. (1) is solved numerically to derive the whole 7r-band structure 
and the coefficients for the TB wave functions. The low-energy band structure around the 
K point in the 1st Brillouin zone is presented in Fig. 2. In this region, the energy dispersion 
can be described by Dirac theory [54]. At zero fields (Fig. 2(a)), all energy bands are spin 
degenerate. The main effect of the SOI is to produce a small energy gap E g between the 
valence and conduction bands. When a perpendicular electric held E z is applied, each 
band exhibits the spin-split subbands due to the breaking of the inversion symmetry (Figs. 
2(b)-2(d)). Consequently, the system is characterized by two energy gaps, Ei g and E 2g 
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for spin-down (blue curve) and spin-up (red curve) subbands, respectively. The former 
decreases with increasing field strength, while the opposite is true for the latter. The 
spin-down and spin-up subbands are reversed at different valleys. Whenever E z = E c 
(~ A so/l — 17 meVA *), the staggered sublattice potential energy exactly compensates 
for the SOI on the spin-down states and the gap of E\ g is closed (Fig. 2(c)). Meanwhile, 
the spin-split energy at the K point reaches its maximum, i.e., 2A so- For E z > E c (Fig. 
2(d)), Ei g is reopened, and the system is transformed into a band insulator. After that, 
both E\ g and E 29 increase with E z , while the spin-split energy at K point remains constant 
at 2 A 50 

The unscreened-excitation spectrum is helpful in understanding the single-particle ex¬ 
citation (SPE) channels. At zero temperature (Fig. 3(a)), the interband excitations are 
the only available excitation channel. The threshold excitation energy is expressed as 
UJ th ter = 2(A so) 2 + Vpq 2 /4:. The real (ti) and imaginary (e 2 ) parts of the longitudinal 
dielectric function, related through the Kramers-Kronig relations, exhibit a logarithmic- 
divergent peak and a step-like structure at respectively. The interband-gap tran¬ 

sitions account for these prominent structures. There are no zero points in e\. That is, 
for T = 0, no self-sustained plasmon excitations can exist in the low-frequency region. 
The intraband transitions are allowed at finite temperatures; these transitions create ad¬ 
ditional special structures in e, i.e., an asymmetric dip and peak in and e 2 , respectively, 
as illustrated in Fig. 3(b) for T=50 K. At the same time, the interband-transition-related 
structures become weaker. The highest intraband-excitation energy is uj™ tra = VfQ, which 
is lower than u;^' ter . Therefore, a SPE gap is developed in e 2 . The asymmetric dip in e\ 
creates two zero points. One of the zero points is located in the e 2 gap (only for the con- 
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dition ksT < E g ), where an undamped plasmon mode can occur. The temperature range 
that allows the undamped plasmon depends on the width of the SPE gap; this width is 
tunable by an external E-held. 

The effects of the E-held on e are displayed in Figs. 3(c)-3(e) for fixed T = 50 K and 
different E z s. At E z = 0.5 E c , shown in Fig. 3(c), the step structure (peak) in 62 (ei) is 
split into two owing to the lifting of the spin degeneracy, as indicated by the pair of blue 
arrows. The two steps of the structure develop at and wgr, respectively. Between 

ujl i'th' an d u Tth T the density of the e-h pairs is half of that beyond ^th- This means 
the creation of a low-damped plasmon region. oj\ n ^ r decreases when E z is increased from 
zero to E c . i.e., the SPE gap is reduced. This may hinder the occurrence of an undamped 
plasmon. The opposite is true for E z > E c (Fig. 3(e)). At E z = E c (Fig. 3(d)), uj™£ r and 
u)™ tra merge. Therefore, both zero points of e\ are located where 62 is finite. Under such 
condition, only damped plasmons can exist. 

The energy-loss function (oc Im[— 1/e]), related to the screened excitation spectrum, 
is used to understand the collective-excitation mode of the electrons. The quantity can 
be measured by the electron energy loss spectroscopy in transmission [50] jEQ 155] [56] and 
reflection m, and inelastic light scattering [5S]. Each prominent structure in Im[—l/e] 
may be viewed as a plasmon excitation with different degrees of Landau damping. At 
T=0, the low-frequency Im[— 1/e] appears as plateaus without and with an electric held, as 
shown in the insets of Figs. 4(a) and 4(b), respectively. The featureless plateau structures 
correspond to the interband excitations, which exhibit single-particle like behaviors. At 
finite temperatures, a prominent peak is induced with a frequency lower than the plateaus 
(Figs. 4(a)-4(d) for T = 50 K). The corresponding dielectric functions reveal that this 


peak could arise from the interplay between the intraband and interband-gap transitions 
(Figs. 3(b), 3(c), and 3(e)) or merely from the sufficiently strong intraband transitions 
under the condition of a zero SPE gap (Fig. 3(d)). The former leads to the stronger peak 
intensity and the narrower peak width in the energy-loss function (Figs. 4(a), 4(b), and 
4(d)), compared with that caused by the latter (Fig. 4(c)). 

The g-dependent behavior of plasmons is important in understanding their character¬ 
istics. Figs. 5(a)-5(c) show the results for different E' z s at fixed T = 50 K. The color scale 
stands for the values of the energy-loss function. The dashed curves in each figure represent 
the highest intraband (white curves) and the lowest interband (blue curves) SPE energies, 
which obey the relations u — VFq and u = 2\J ( £E Z — s^/X | G + a 2 X 2 R q 2 /4 ) 2 + h 2 Vpq 2 /4 
with s = ±1, respectively. The expressions are derived from the analytical solution of 
single-particle spectrum near the K point [M]. These boundary curves define the various 
regions with different degrees of Landau damping. Between the white and the blue bound¬ 
ary curves lies the SPE gap, where exists an undamped plasmon branch. The plasmon 
frequency u p is well fitted by yfq in the long-wavelength limit, a feature of a 2D plasmon. 
For lower temperatures (Figs. 5(a) and 5(b)), this plasmon mode is finally damped out in 
the intraband SPE continuum at larger q's. The existence range of the undamped plas¬ 
mon in momentum is proportional to the width of the SPE gap, which is reduced in the 
region of 0 < E z < E c (Fig. 5(b) for E z = 0.5 E c ). There are two blue boundary curves 
in E z 7 ^ 0 (Figs. 5(b)-5(d)). They are the result of the lifting of the spin degeneracy. 
Between the two blue curves lies a low damped region. In Fig. 5(b), a broad plasmon 
branch appears in the region, just above the undamped plasmon mode. The coexistence 
of the two collective-excitation modes can only be realized at ksT m E Xg . The SPE gap 
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is closed at E z = E c (Fig. 5(c)). Then, the plasmon is damped at arbitrary q's. Until 
E z > E c , the undamped plasmon branch is regained (not shown). Under the above con¬ 
ditions, the plasmon dispersion is confined between two SPE boundary curves. However, 
when the temperature is sufficiently high, e.g., T = 200 K and E z = 0.5 E c in Fig. 5(d), 
the plasmon branch crosses the SPE boundaries and shows abrupt changes in intensity 
and bandwidth. The discontinuous plasmon dispersion is more evident when the plasmon 
crosses the upper blue curve. After that, the plasmon enters the strong damped region. It 
is worth noting that the temperature-induced plasmon differs from the extrinsic plasmon 
discussed in doped silicene |5PJ. The dispersion relation of the extrinsic plasmon is decided 
by the Fermi energy. This is a great contrast to the intrinsic plasmon whose dispersion is 
dominated by the temperature. The different electron distributions result in the different 
plasmon dispersions between in the two conditions. 

The effect of temperature variation on the plasmon excitation is worth a further discus¬ 
sion. Figs. 6(a)-6(c) show the T-dependent plasmon spectra at fixed q = 1 and different 
E' z s. The distribution range of the undamped plasmon mode (between the blue and white 
dashed curves) in T is reduced in the region of 0 < E z < E c (Fig. 6(b) for E z = 0.5 E c ) 
and diminished to zero at E z = E c (Fig. 6(c)), compared to the condition of E z = 0 (Fig. 
6 (a)). The undamped plasmon is transformed to a damped one above a critical tempera¬ 
ture (T c ~ 100 K for E z = 0). The transformation in the plasmon mode during the change 
in T is a discontinuous process, a feature that could allow for the identification of the spin- 
orbit energy gap. The discontinuous dispersions of the plasmons remain when changing 
background dielectric constant eo and energy broadening T. In Fig. 6(b) for E z = 0.5 E C) 
two critical temperatures exist through which the plasmon spectrum is drastically altered. 
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The two T' c s correspond to the two energy gaps opening in the system. This conlcl be 
used to examine the spin degeneracy and determine E c . Note that above the highest T c 
(~ 175 K for E z = 0.5 E c ), the plasmon intensity grows monotonically with T, similar to 
the behavior of the intrinsic plasmon in gapless graphene [T3]. This indicates that at such 
high temperatures the SOI effect becomes negligible. 

For fixed q and T, the plasmon existence and its spectral weight can be effectively 
controlled by varying E z , as shown in Fig. 7. At T = 50 K and q — 1 (Fig. 7(a)), the 
undamped plasmon is stable for E z < 0.5 E c due to the significantly large SPE gap. When 
E z is increased to near E c , the plasmon is damped by the interband SPEs and faints. 
The SPE gap is reopened in the region of E z > E c . After that the plasmon is intensified 
first with the increment of E z , and when E z exceeds another critical value E' c (E' c « 3 E c 
for T = 50 K and q — 1), the plasmon is gradually damped out in the intraband SPE 
continuum. At T = 100 K (Fig. 7(b)), the undamped plasmon barely exists in E z < 2 E c , 
because the overly strong intraband transitions make the plasmon frequency higher than the 
threshold interband SPE energy. Until E z > 2 E c , the energy gap becomes larger than fcgT, 
and the undamped plasmon emerges. The undamped plasmon at T = 100 K can survive 
in a larger E z than that at T = 50 K. The stronger intraband transitions can explain 
this. At T = 200 K (Fig. 7(c)), the plasmon spectrum displays three distinct degrees of 
Landau damping in different E z ranges. In the first region (0 < E z < 0.7 E c ), the plasmon 
is located beyond the two interband SPE boirndaries and suffers strong Landau damping. 
In the second region (0.7 E c < E z < 2.5 E c ), the plasmon lies between the two interband 
boundaries and is weakly damped. In the third region (E z > 2.5 E c ), the plasmon lies inside 
the SPE gap and is undamped. Though the plasmon intensity strongly depends on E z , the 
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change in the plasmon frequency is insignificant. This illustrates that the appearance of 
the low-frequency plasmon is mainly dominated by the intraband transitions. 

IV. SUMMARY AND CONCLUSIONS 

The Coulomb-excitation properties of silicene are studied by the tight-binding model 
and random-phase approximation. The temperature and E-field effects are considered 
simultaneously. At finite temperatures, an undamped plasmon is induced by the interplay 
between the intraband and interband-gap transitions. The plasmon mode only exists in 
a limited region of T and q, which can be further tuned by varying E z . Thanks to the 
opening of two energy gaps from the spin split in E z , there are three distinct degrees of 
Landau damping in different q, T, and E z ranges. Therefore, the plasmon dispersion is 
discontinuous at certain critical q's and T's. The phenomenon could be used to identify the 
spin degeneracy and to determine relevant theoretical parameters, such as the magnitude 
of SOI. We further show that, for given q and T, changing E z can control the plasmon 
intensity and its existence at a nearly fixed frequency, a property that should be important 
in plasmonic applications. The tunable low-frequency plasmons can potentially be also 
found in other 2D crystals with band gaps that can be adjusted by an external E-field, e.g., 
germanene and stanene. 
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Figure Captions 


Figure 1. The lattice geometry of low-buckled silicene. The two sublattices are sepa¬ 
rated by a perpendicular distance 21. The staggered sublattice potential energy is produced 
by an external electric field E z . 

Figure 2. The low-energy bands of monolayer silicene in the presence of spin-orbit 
coupling about the K point are shown in (a). Those under E z = 0.5 E c , E c , and 2 E c are 
shown in (b), (c), and (d), respectively. 

Figure 3. The dielectric function as a function of oj at q = 1 and T = 0 is shown in (a). 
Those for T = 50 K and different E' z s are shown in (b)-(e). The purple and orange curves 
represent the real and the imaginary parts, respectively. The energy broadening width is 
T = 0.1 meV. q is in unit of 10 5 /cm. 

Figure 4. The energy-loss functions for fixed q — 1, T — 50 K, and different E' z s are 
plotted in (a)-(d). Those for E z — 0 and E z = 0.5 E c at T = 0 are plotted in the insets of 
(a) and (b), respectively. 

Figure 5. The q-dependent plasmon spectrum at T = 50 K for (a) E z = 0, (b) E z = 
0.5E C , (c) E z = E C1 and (d) E z = 2 E c . The dashed blue and white curves in each panel 
indicate the interband and intraband SPE boundaries, respectively. The color scale stands 
for the values of the energy-loss function. 

Figure 6. The T-dependent plasmon spectrum at fixed q — 1 and different E z s: (a) 
E z = 0, (b) E z = 0.5 E c , and (c) E z = E c . The blue and white dashed curves in each panel 
indicate the interband and the intraband SPE boundaries, respectively. The color scale 
stands for the values of the energy-loss function. 
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Figure 7. The independent plasmon spectrum at fixed q — 1 and different T's: (a) 
T = 50 K, (b) T = 100 K, and (c) T = 200 K. The color scale stands for the values of the 
energy-loss function. 
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